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The HI Collaboration recently reported a new analysis on the inelastic photo-production of J /ip 
mesons at DESY HERA ep collider. We show that these new experimental results are well described 
by the Color Evaporation Model for quarkonium production. Moreover, this new data requires the 
introduction of resolved photon contributions in order to explain the results on small charmonium 
energy fraction, indicating that colored cc pairs also contribute to the process. 
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I. INTRODUCTION 


The HI Collaboration recently reported an analysis on 
the inelastic photo-production of J/if> mesons jp where 
their new data was confronted with the color singlet j|] 
and color octet || models. We compare the HI data 
with the color evaporation model (CEM) predictions for 
quarkonium production showing that it provides a good 
description of this new set of data. Moreover, this new 
data probes for the first time the small charmonium en¬ 
ergy fraction (z) region, requiring the introduction of re¬ 
solved photon contributions in order to explain the re¬ 
sults. These contributions are basically due to color octet 
configurations, confirming the correctness of their inclu¬ 
sion in the charmonium production mechanism. 

Previous measurements of the inelastic photo¬ 
production of charmonium at HERA 1 appeared 
to have ignited a charmonium crisis. The color singlet 
model (CSM) to onium production cross section fits the 
data for large charmonium energy fraction z, where color 
octet models seemed to fail. This fact is, however, in 
qualitative disagreement with a wealth of information 
that exists on charmonium production by other initial 
states. We suggested [(J that this discrepancy is due 
to the neglect of non-perturbative effects that are im¬ 
portant at large z. Implementing a phenomenological 
parametrization of these effects in a scheme originally 
developed for Drell-Yan phenomenology, we illustrated 
how agreement with data could be achieved. In this 
work, we employ the procedure of 0 to mimic the non- 
perturbative effects. 


II. THEORETICAL FRAMEWORK 


It is clear nowadays that J/•ip production is a two-step 
process where a heavy quark pair is produced first, fol¬ 
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lowed by the non-perturbative formation of the colorless 
asymptotic state. The predictions of the once conven¬ 
tional treatment of color in perturbative QCD calcula¬ 
tions, i.e., the CSM 0], for the charmonium production 
at the Fermilab Tevatron is at variance with the avail¬ 
able data |0| . As a consequence, color octet as well as sin¬ 
glet cc states contribute to the production of J/ip. Two 
formalisms have been proposed to incorporate these fea¬ 
tures: the Non-Relativistic QCD, also known as COM 
||, and the renewed CEM scheme M [TT) , also known 
as Soft Color Interactions (SCI) jid|. The original CEM 
||T0 actually predates the color singlet approach, and had 
been abandoned for no good reason. Recent measure¬ 
ments of the polarization of bound charm and beauty 
mesons, seems to disfavor the COM framework fll4||. 


The color evaporation model simply states that char¬ 
monium production is described by the same dynamics 
as DD production, i.e., by the formation of a colored 
cc pair. Rather than imposing that the cc pair is in a 
color-singlet state in the short-distance perturbative di¬ 
agrams, it is argued that the appearance of color-singlet 
asymptotic states solely depends on the outcome of large- 
distance fluctuations of quarks and gluons. These large- 
distance fluctuations are probably complex enough for 
the occupation of different color states to approximately 
respect statistical counting. In other words, the for¬ 
mation of color-singlet states is a non-perturbative phe¬ 
nomenon. In fact, it does not seem logical to enforce the 
color-singlet property of the cc pair at short distances, 
given that there is an infinite time for soft gluons to read¬ 
just the color of the pair before it appears as an asymp¬ 
totic if}, Xc or, alternatively, DD state. It is indeed hard 
to imagine that a color-singlet state formed at a range 
nr))) 1 automatically survives to form a if). This approach 
to color is also used to formulate a successful prescrip¬ 
tion for t he p roduction of rapidity gaps between jets at 


Tevatron |7|, |l8|, |l9j, |2Cj, |2l|] and HERA 


Although far more restrictive than other proposals, the 
CEM successfully accommodates all features of charmo¬ 
nium production 100- It predicts that the sum of 
the cross section of all onium and open charm states is 
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where M cg is the invariant mass of the cc pair. The factor 
1 /9 stands for the probability that a pair of charm quarks 
formed at a typical time scale 1 /M$ ends up as a color 
singlet state after exchanging an uncountable number of 
soft gluons with the reaction remnants. One attractive 
feature of this model is the relation between the produc¬ 
tion of charmonium and open charm which allows us to 
use the open charm data to normalize the perturbative 
QCD calculation, and consequently, to make more accu¬ 
rate predictions for charmonium cross sections. 

The fraction p ^ of produced onium states that mate¬ 
rialize as ip, 


G'lfj — P-0 ^onium 5 


( 3 ) 


has been inferred from low energy measurements to 
be a constant (24|, [|^]. From the charmonium photo¬ 
production, we determined that p^, = 0.43-0.5 Q; a 
value that can be accounted for by statistical counting 
of final states [j^. The fact that all ip production data 
are described in terms of this single parameter, fixed by 
J/ip photo-production, leads to parameter free predic- 
tions for Z-boson decay rate into prompt ip | 2 (| , and 
to charmonium production cross section at Tevatron |27j] 
and HERA ||, as well as in neutrino initiated reactions 
p 8 [ . These predictions are in agreement with the avail¬ 
able data. 


III. INELASTIC PHOTO-PRODUCTION OF 
CHARMONIUM 


Application of the CEM scheme to inelastic charmo¬ 
nium photo-production is straightforward. According to 
the parton model, the cross section for the J/ip photo¬ 
production at a given center-of-mass energy W is 


^ 7 p — >J/ipx(A ) II pA/'y( x A ) fB/p( x B ) dAB — ^J/i/jX( s') dXA dXB 

I- 


( 4 ) 


where the subprocess cross section a is given by Eqs. (|l]) 
and ©• Here, = y/ xaXb W is the center-of-mass 
energy of the subprocess AB —> J/ipX , and /a / 7 (/b/ p ) 
is the distribution function of the parton A ( B ) in the 
photon (proton). For direct photon interactions (A = 7 ) 
we have /a/^xa) = d(x A - !)• 

An important kinematical variable is 


0 = 


Pj/P ‘ Pp 

p . p 

1 p 


( 5 ) 


where Pj/p t -y tP is the four-momentum of the J/ip, photon, 
or proton, respectively. In the proton rest frame, z is the 
fraction of photon energy carried by the J/ip. 

The lowest order CEM contribution for the charmo¬ 
nium production is the direct photon process 7 g —> cc. 
However, it is important only for z ~ 1. For the range 
of z we are interested, the direct photon contribution is 
dominated by the diagrams depicted in the Fig. 0(a)- 
The charm quark pair in 7 g fusion can be produced 
in both color singlet and octet configurations, while 7 q 
fusion leads only to colored cc pairs. Besides the di¬ 
rect photon-gluon and photon-quark reactions, we also 
included resolved photon processes, which proceed via 
quark-quark, quark -gluon, and gluon-gluon fusion into 
cc+quark (gluon), as exemplified in Fig. |I](b). 
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FIG. 1: Processes contributing to direct (a) and resolved (b) 
charmonium photo-production. 


The higher order processes, like the ones in Fig. |l], have 
to be evaluated with some caution in the region of rela- 
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tively small J/?/> transverse momentum and large z. For 
small f-channel momentum transfer (t —> 0 ), the gluon 
exchange diagrams in Fig. [i] represent the QCD evolution 
of the initial state gluon distribution functions, and not 
higher order gluon exchange diagrams. Therefore, the 
contribution of the gluon-exchange diagrams where t is 
soft, has already been taken into account by the leading 
order 7 g —> c c diagram |2tJ. By itself, when evaluated 
at tree level, the diagrams in Fig. lead to a divergence 
for z —> 1 , that appears as an unphysical growth of the 
cross section for z < 1. Although the double counting is 
relevant only for vanishing pp and z = 1 , the enhance¬ 
ment of the cross section associated with the collinear 
divergence will, in practice, affect the calculation until 
Pt and z-values become truly perturbative. 

Presently, complete higher order QCD calculations 
are not available for the processes we are interested 
in. Therefore, we phenomenologically took into ac¬ 
count non-perturbative QCD corrections in large z char- 
monium photo-production by introducing the following 
parametrization to mimic the non-perturbative contri¬ 
butions |ij : 

5 gj=[l-FW». W )][l-G«3„, 2 )]gg (6) 

with 

_** 

F(Q 0 ,p T )=e ^ (7) 

and 

G(Q 0 ,z) = e~^, (8) 

where atree is the tree level perturbative cross section, 

= z(Ql+4ml) — Aml, and (1 — z 0 ) = {p^+4m 2 )/ (Q§-F 
4m 2 ) are positive definite or null. The parameter Q 0 in¬ 
dicates the value of the momentum transfer t where the 
higher order diagrams contribute in the truly perturba¬ 
tive regime |3(J. We found in |b| that Q 0 = 2 m c best 
describes the data, and it is this value that we used in 
the present work. 


IV. RESULTS 


In this work we evaluated numerically the tree level 
scattering amplitudes using the Madgraph 131 and Helas 
fi2| | packages. The phase space integration was per¬ 
formed using the adaptative Monte Carlo program VE¬ 
GAS ia. 


There is a wide choice of QCD parameters appearing in 
the evaluation of charmonium photo-production, which 
leads to a large theoretical uncertainty. For the sake of 
definiteness, we present only the uncertainty that arises 
from possible choices of the charm quark mass, keeping 
all others QCD parameters fixed. All our results are pre¬ 
sented as a central value surrounded by an error band 
corresponding to the choice of the charm quark mass as 


to c = 1.3±0.1 GeV. We used the GRV-94 LO Q param¬ 
eterization of the proton structure functions and GRV-G 
LO Q for the photon parton density. For both structure 
functions, we set the factorization scale as pf = V§- We 
evaluated the running of the strong coupling constant in 
leading order with four active flavors using A qcd = 300 
MeV, renormalization scale pr = \/2m c for direct pro¬ 
cesses, and pr = V§ for resolved processes. We also used 
= 0.5 as the J/ip fraction of all charm bound states. 

The HI Collaboration performed their analysis using 
data in which Q 2 < 1 GeV 2 , and subdivided their data 
into several different kinematical regions, in order to bet¬ 
ter determine the region where perturbative QCD calcu¬ 
lations furnish a reliable description of the data. 




W yp (GeV) 


FIG. 2: Total cross section as function of W lp for pr > 1 
GeV: (a) 0.3 < 2 < 0.9 and (b) 0.3 < z < 0.8. The shaded 
band shows the theoretical prediction for m c = 1.3±0.1 GeV. 


We collated the CEM predictions with the HI results 
using the cuts applied to the experimental data. Initially 
we analyzed the behavior of the total cross section as 
function of the "fp center-of-mass energy (W) for two 
different z regions, requiring a minimum J/ij) transverse 
momentum pr > 1 GeV. In Fig. ||(a) we show the CEM 
predictions for 0.3 < z < 0.9, and in Fig. |(b) we present 
our results for 0.3 < 2 < 0.8. As we can see from these 
figures, moving the upper limit in z alone does not change 
much the shape of the curve, and within theoretical un¬ 
certainties, agreement is found for both regions. 

We display in Fig. ||, the CEM predictions for the z 
spectrum, requiring pr > 1 GeV and a reaction center 
of-mass energy in the range 120 < W lv < 260 GeV. In 
this figure, the experimental points represented by trian- 

B 1 3 S and squares stand for two different datasets; see Ref. 

for details. The direct photon contribution is repre¬ 
sented by the dashed line, and the resolved one by the 
dotted line; the solid line displays the sum of direct and 
resolved contributions. The shaded area corresponds to 
the total CEM prediction for the charm quark mass de¬ 
scribed above. As we can see from this figure, the agree¬ 
ment of the CEM with data is quite good. Moreover, the 
data at low and medium z do require the introduction of 
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FIG. 3: Differential cross section as function of the inelasticity 
parameter 2 for p^ > 1 GeV 2 and 120 < W lp < 260 GeV. The 
dashed line stands for the contribution from direct processes 
while the dotted line shows the contribution from resolved 
processes. The solid line is the sum of both contributions. 
The shaded band shows the theoretical uncertainty on the 
value of the charm mass (m c = 1.3 ± 0.1 GeV). 

resolved processes in order to explain the results. This is 
a clear signal that colored charm quark pairs contribute 
to the J/if.’ production. 

In Fig. [| we present the p\ distribution for the low z 
sample 0.05 < z < 0.45 which corresponds to the trian¬ 
gles in Fig. |j. For these values of z we expect the theo¬ 
retical uncertainties due to higher order corrections to be 
small, and we can see that the CEM and data agree well. 
Consequently, the disagreement between data and theory 
at large z is indeed due to the importance of higher order 
non-perturbative corrections. 



FIG. 4: Squared transverse momentum p^ distribution for 
data collected in the very inelastic region 0.05 < z < 0.45 
and 120 < W lv < 260 GeV. The shaded band shows the 
theoretical uncertainty on the value of the charm mass (m c = 
1.3 ±0.1 GeV). 

Higher order effects are sizeable on the medium and 
high z regions. Therefore, the HI Collaboration divided 
their data sample into several pp and z regions for a 


center-of-mass energies in the range 60 < W lp < 240 
GeV. This allows us to better compare data with theory. 



z 


FIG. 5: Differential cross section as function of the inelas¬ 
ticity parameter z for pp > 1, 2, and 3 GeV, divided by a 
factor 1, 10, and 100 respectively. The shaded bands show 
the theoretical prediction obtaining varying the charm mass 
(m c = 1.3 ± 0.1 GeV). 

In Fig. we compare the CEM predictions for the z 
spectrum with data for three different values of the J/ip 
minimum transverse momentum, i.e., pp > 1, 2, and 3 
GeV. The curves have been divided by factors 1, 10, and 
100, respectively, in order to help visualization. There 
is an overall agreement between data and theory, ex¬ 
cept for the highest bin in z. Moreover, removing our 
regularization procedure only worsens the theoretical re¬ 
sults [Q . Increasing the minimum value of the transverse 
momentum does not improve the quality of the fitting, 
what means that our parametrization of the higher or¬ 
der effects has correctly incorporated the dependence on 
minimum pp. 



FIG. 6: Squared transverse momentum p P distribution for 
0.3 < z < 0.9 and center-of-mass energy 60 < W lp < 240 
GeV. The shaded band shows the theoretical uncertainty on 
the value of the charm mass (m c = 1.3 ± 0.1 GeV). 

Fig. 0 contains the J/ip squared transverse momentum 
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FIG. 7: Squared transverse momentum p\ spectrum for 60 < 
W-yp < 240 and data collected in three different inelastic bins: 
0.75 < a < 0.9 (upper), 0.6 < a < 0.75 (middle), and 0.3 < 
2 < 0.6 (lower). To help visualization the curves have been 
scaled as 1, 1/10, and 1/100 from top to bottom. The shaded 
bands show the theoretical uncertainty on the value of the 
charm mass (m c = 1.3 ± 0.1 GeV). 

distribution taking into account the medium and high z 
data. The agreement between theory and data is quite 
satisfactory. However, the shape of the spectrum shows 
some disagreement at the very low pp bins. In order to 
understand what is happening, let us consider this distri¬ 
bution for three different z bins, namely 0.75 < z < 0.9, 
0.6 < 2 < 0.75, and 0.3 < 2 < 0.6; see Fig. || where the 
curves have also been divided by factors 1 , 10, and 100, 
respectively, to help visualization. We can learn from this 
last figure that the agreement between CEM predictions 


and data improves for low 2 regions. This fact can be 
understood as a limitation of the proposed parameteriza¬ 
tion to completely mimic higher order QCD contributions 
when we approach the elastic region. It also implies that 
the observed small discrepancy with data for the highest 
bin on 2 must be credited to the lack of a complete QCD 
calculation and is not related to the Color Evaporation 
approach to describe quarkonium production. 


V. CONCLUSIONS 

We showed that the Color Evaporation Model de¬ 
scribes the available data on J/i/> photo-production, pro¬ 
vide we include high order QCD corrections at high in¬ 
elasticities 2 . Moreover, the newly available data at low 
2 provides a clear proof of the importance of colored cc 
pairs to the production of charmonium, since the data 
on this region can only be explained considering resolved 
photon processes, which lead to colored cc pairs. 
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